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Abstract 

By means of the Brownian dynamics (BD) method of simulations we have developed, we study 
dynamics of individual DNA undergoing constant field gel electrophoresis (CFGE), focusing on 
the relevance of the 'defect' concept due to de Gennes in CFGE. The corresponding embodiment, 
which we call slack beads (s-beads) is explicitly introduced in our BD model. In equilibrium under 
a vanishing field the distance between s-beads and their hopping range are found to be randomly 
distributed following a Poisson distribution. In the strong field range, where a chain undergoes the 
elongation-contraction motion, s-beads are observed to be alternatively annihilated in elongation 
and created in contraction of the chain. On the other hand, the distribution of hopping range 
of s-beads does not differ much from that in equilibrium. The results indicate that the motion 
of the chain elongated consists of a huge number of random movements of s-beads. We have 
also confirmed that these features of s-beads agree qualitatively with those of s-monomers in the 
extended bond fluctuation model (EBFM) which we recently proposed. The coincidence of the 
two simulations strongly supports the stochastic semi-local movement of s-monomers which we a 
priori introduced into the EBFM. 
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I. INTRODUCTION 



Gel electrophoresis is a major technique which separates polyelectrolytes according to 
their size. Among them gel electrophoresis of DNA is one of the most interesting problems 
of statistical physics of polymers. Computer simulationsBIl and experimentsl'S reported 
a peculiar behavior of individual DNA that they exhibit contracted and extended forms 
alternatively under a constant field. Recently we have developed a Brownian dynamics 
(BD) method which simulates electrophoresis of DNA in a 3 dimensional (3D) space by a 
chain of electrolyte beads. Its details have been reported in our preceding paper (hereafter 
refereed to as I), where we have demonstrated 'quasi-periodic' time evolution of the velocity 
of the center of mass, VG{t), and of the radius of the longer principal axis, Ri{t), of the chain. 
The present paper is devoted to clarify the role of 'defects' introduced by de Gennesi on the 
elongation-contraction motion. 

The concept of 'defect' proposed by de Gennedl is useful to investigate dynamics of 
individual polymers confined to a space with entanglement such as gel electrophoresis of 
DNA. A defect represents a relatively loose portion, or a slack segment, of a polymer and 
is assumed to behave as a particle of ideal gas. He showed that its diffusion in the tube 
made of other polymers leads to the famous growth law 0(t) oc t^^^ in the short time regime, 
where (f){t) is the mean square displacement of a segment in equilibrium. 

Electrophoresis of a chain in an ideal regular network of other fixed chains has been 
extensively studied by Monte Carlo (MC) simulations.S0'i Among them, the repton model 
proposed by Rubinstein is important.! It simulates dynamics of a chain only by stochastic 
movements of reptons each of which can be regarded as embodiment of a defect. The 
model is applied to simulation of the constant field gel electrophoresis (CFGE) by Duke.0 
It was shown that the mobility is inversely proportional to the number of reptons in a very 
weak field region. The relation was consistent with what the reptation theory predicted 
under the vanishing field limit. Under finite fields, however, the results derived from 
the repton model qualitatively differ from those derived from the biased reptation theorylll: 
the mobility plateau is reached at lower fields in the former model than those predicted by 
the latter theory. Duke also observed snapshots of tubes which look like hooks, and argued 
that once a chain is hooked around a gel fiber, both its arms tend to stretch and rarely 
retract. This is because defects prefer to move in direction of the field and to create new 
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tube segments which in turn prefer to ahgn to the field. The observed fact that the chain 
rarely exhibits a retraction is considered to be an artifact of the MC method which uses 
local hopping rules. The incorporation of a hopping rule which allows non-local movement 
of defects was done by Duke and Viovy in their advanced version of the repton model.lii0 
Simulations of chains with over a thousand segments revealed that a chain exhibits ramified 
configurations where 'hernias' are created and annihilated under a constant field.lii 

Recently, we have investigated gel electrophoresis and diffusive dynamics of a chain in a 
space with immobile obstacles by making use of an extended version of the bond fiuctuation 
model (BFM).El The conventional BFM (CBFM)0 is known as an efficient model to simulate 
general problems involved in polymers. However, it is faced with a problem which is similar 
to that of the repton model mentioned above: the chain motion extremely slows down 
when it is applied to CFGE under relatively large fields.! To overcome this problem, we 
have empirically incorporated semi-local dynamics of slack parts of a chain. The latter 
correspond to 'defects' due to de Gennes and were named s- monomers. They are explicitly 
defined in each configuration appearing in the CBFM. We have called such a BFM with 
introduction of semi- local movements of s- monomers the extended BFM (EBFM). Both the 
DV model and the EBFM incorporate nonlocal movements of slack portions of DNA. There 
are, however, the following difference between the two methods. In the DV model any defect 
is not created or annihilated except at the extremity, whereas in the EBFM, s-monomers 
are continually created and annihilated anywhere in the chain, since it involves the CBFM 
trials of conventional local movements of all monomers. 

Using the EBFM it has been confirmed that the trapping problem of the CBFM under 
large fields is overcome and the field-dependence of mobility is reproduced qualitatively as 
observed by the experiment.0 Furthermore, quasi-periodic behavior, or a damped oscilla- 
tion, is observed in the correlation function of Ri{t) in 3D0 as it has been observed in the 
experimentiS as well as by the BD calculations in 2D (without the excluded volume effect 
and in 3D (with the effect). 

The defect is thus a quite important concept to understand DNA dynamics in gel. The 
purpose of the present paper is to clarify the microscopic bases of the defect dynamics, 
thereby presenting a microscopic support of the EBFM. For this purpose we examine the 
corresponding object in a chain, i.e., slack beads (abbreviated as s-beads) by means of the 
BD method. The work which is similar to this approach has been done by Masubuchi et al.iil 
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in their 2D BD study on the biased sinusoidal field gel electrophoresis (BSFGE). Looking at 
time evolution of snapshots of kinks which correspond to slack segments, they have found 
that a portion with less kinks develops near the segment which is hindered by an obstacle, 
and that the portion gradually shifts as the chain slides off the obstacle. They have argued 
that the growth of kinks plays a key role on the decrease of the mobility observed in chains 
with particular sizes under BSFGE. 

In the present paper we examine static and dynamic properties of defects much more 
in details by means of our BD method for a chain in a 3D space. As described in our 
accompanying paper, hereafter referred to as 1,0 the method incorporates semi-microscopic 
ingredients of dynamics of a polymer, i.e., the excluded volume effect, the non-linear elastic 
force which keeps the distance between the two adjacent beads (segments) and the random 
force from a solvent. In each configuration of the chain simulated by the BD method s- 
beads, or 'defects', are explicitly specified in a quite similar way as the s-monomers have 
been defined in the EBFM. Then we investigate various statistics associated with the s- 
beads, or more explicitly, portions between the neighboring s-beads which we call extended 
parts of the chain. 

In particular, we have concentrated on the normalized histogram of the extended parts 
of length n, rn{t), and that of their displacement along the chain direction, rAu(t). In 
equilibrium under a vanishing field it is found that both r„(t) and r/^u{t) obey Poisson dis- 
tributions as expected. In the elongation-contraction motion under CFGE, r„(t) exhibits 
a peculiar time evolution refiecting the elongation and contraction of the chain. The his- 
togram r/\uit), on the other hand, does not fiuctuate correspondingly to r„(t) but exhibits 
random hopping nature rather similar to the one observed in equilibrium. These results 
indicate that defects move quite stochastically even in the deterministic regime where the 
chain elongation-contraction motion looks apparently deterministic. These features of r„ (t) 
and rAu(t) are confirmed to be in qualitative agreement with those observed for s-monomers 
in the EBFM. 

The organization of the present paper is as follows. In the next section we briefiy explain 
our numerical method, and discuss equilibrium dynamics of s-beads in Sec. |T|. Charac- 
teristic features of r„{t) and r^u{t) in the field range where a chain exhibits elongation- 
contraction motions are examined in Sec. |V|. The corresponding results obtained by the 
EBFM are discussed in Sec. |V|, and the final section is devoted to the conclusion. 
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II. NUMERICAL METHOD 



The BD method we adopt in the present work is to solve the following Langevin equation 
of motion for a chain of spherical beads in a continuous 3D space, 

a j a 

where Xi is the position of the i-th bead, and (, q, and Ey, are the viscosity, the charge of 
a bead and the (bare) external field, respectively. The second and the third terms in r.h.s. 
are the constraining forces required for the beads to satisfy the conditions, 

||Zi|| = \\xi+i - Xi\\ < \p2. (3) 

putting the bead diameter unity. The last term in eqn. ^ stands for the random force from 
a solvent. Explicitly we adopt the following form for it 

m) = Y,n8{t-n^t), (4) 

n 

where At is regarded as the average period of random forces. The distribution of is 
characterized by the averages 

(/D = (5) 
(/r/D = 2CfcBTAt<5,,5^„. (6) 

Here ks and T are the Boltzmann constant and the absolute temperature, respectively. 

The details of solution of the above set of equations have been explained in our accom- 
panying paper I. As a representative of 'defect', or a loose portion of a chain in the BD 
calculation we define a slack bead (s-bead). For any sequences of three beads, if the dis- 
tance between the first and the third is not greater than v^, the second one is regarded as 
an s-bead. This definition of 'defect' is the same as that of s-monomers in a 3D version of 
the EBFM we investigated before@ (see in Sec. |Vp. We call a sequence of beads between 
neighboring s-beads an extended part of the chain. 

III. 'DEFECTS' IN EQUILIBRIUM 

In equilibrium under a vanishing field, 'defects' are considered to be distributed randomly 
along the polymer coordinate. Correspondingly, in our BD analysis, s-beads are randomly 



distributed throughout the chain. Then the 'extended' parts consisting of n beads is expected 
to obey a distribution 

(M-n)(l-p)" 



E£i(M-z)(1-p)^ 



(7) 



where M is the total number of beads of the chain and p the density of s-beads. For 
M ^ 1, reduces to a Poisson distribution: 
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Next let us consider A-u, displacement of the tangential component, of extended parts 
between neighboring s-beads. In our BD calculation, it is evaluated as 

im + l 
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X 



X 



i+1 



(9) 



where Axi is the displacement of the i-th bead in one MD step solving the Langevin equa- 
tion III and index denotes a figure which is numbered on the s-beads. From eqn. |l| we 
may assume that the distribution of Axi is essentially given by that of the random force, 
defined by eqn. M. Then the distribution of Asi is written as 
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■exp(-CAs2/2A;BT). 
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Using this S(Asj) and \&(?t,) of eqn. H, we can evaluate the distribution of Am as 
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Thus <I>(Am) is expected to obey also a Poisson distribution with the width f^An 



1 / ^\og{-^) / [k-oT /2C,). The latter is proportional to y/k^T if the T-dependence of p is 
neglected. 

In Fig. m, we show the normalized histograms r„ and in equilibrium for various T 
which are obtained by the BD calculation. Here r„ = Nn/ Yin where Nn is the number 
that the extended parts of a length n appear in the chain configurations simulated, and tam 
is similarly defined. In agreement with \I^(n) of eqn. ^ r„ obeys a Poisson distribution. The 



exponential tail of r„ at k^T = 1 is fitted to e~"/'^ with u = {n) = 5.7, which yields p = 0.17 
through the relation to = l/log^Y3^)- simulation we have also directly counted 

the number of s-beads, which yields the same s-bead density p = 0.17 as expected. The 
distribution r^u of the tangential movement of the extended parts, which we regard as that of 
the s-beads, also obeys a Poisson distribution as seen in the same figure. Its exponential tail, 
specified by ^Iau as e""''*^"^" , significantly depends on T. As shown in Fig. 0, its dependence 
on T is consistent with eqn. |1^. Even the pre-factor 0.30 of the proportionality to y/ksT 
agrees with the predicted value 1 / ^2 log = 0.37 within 20%. 

The results obtained above confirm the role of the 'defects', or the s-beads, on the static 
and dynamic properties of a polymer in equilibrium. In particular, it is natural to regard 
c{n){l) as the 'persistence length' of the chain, where (n) is the mean value of n and c is 
a numerical constant nearly equal to unity. In this context, we note \/2 > (Z) > 1, where 
(/) is the mean distance between two adjacent beads. In fact we have already argued in I 
that this interpretation yields an appropriate conversion factor between the length scales of 
the BD simulation and the experiment. Here we have to emphasize that the value of (n) is 
almost independent of T in the range k^T < 8, as well as of M and a^ci we have examined. 
It is observed, however, to reduce about by a half when the excluded volume effect between 
beads is discarded. This enlargement of the 'persistence length' is one of the important 
consequences of the excluded volume effect on the present BD chain in a 3D space. 



IV. 'DEFECTS' DYNAMICS UNDER CFGE 

Under relatively strong fields DNA undergo the quasi-periodic behaviori: real time im- 
ages of individual DNA exhibit elongated and contracted shapes alternatively. Correspond- 
ing behavior was observed in the previous 2D BD mo delUJii likewise by our present BD 
simulation which has been discussed in details in I. We show in Fig. 0, a normalized curve 
of Ri{t), i.e., Ri(t — t^ax) / Riituiax) of chains with M = 160 in a lattice of gel with the 
lattice distance of a^ei = 20 under a field E = qEy, = 0.032. In order to get further in- 
sights into the mean elongation-contraction motion through the defect dynamics, we have 
examined histograms of length of the extended parts Nn{t) and of their tangential move- 
ment N^uit) at four representative time intervals, i.e., [tram — 0.1-D,tmm]5 [^min^^mm + O.ID], 



- 0.1L',tmax], and [tmax, W + O.ID], as indicated in Fig. |. Here t^m = ^min:i - 
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and D = f^i^.-^ — tmini' where t"^^^ is the time of local maximum in the n-th peak in Ri{t), 
^min:r ^min:i respectively local minima just before and after the local maximum at 
^'^d the overline stands for the average over peaks in Ri{t). As for the procedure how to 
pick up the peaks in Ri{t), the reader may refer to I. We denote the four intervals introduced 
here as Interval 1, 2, 3 and 4, respectively. It is expected that Interval 1 corresponds to the 
time range where the chain is in a coiled shape, Interval 2 to that of a contracted shape 
when a U-shape conformation has just started to grow due to trapping by gel, Interval 3 to 
that of an elongated U-shape when the chain is going to get rid of the trap, and Interval 4 to 
that of an I-shape when it is going to retrieve another coiled shape. We cut out sequential 
2000 MD steps (mds) in each of the four intervals of each peak in Ri{t) and accumulate data 
associated with them over 80 peaks whose D are about 1.5 x lO^mds or larger. 

The normalized histograms r„(t) (= Nnit) / Y^^Nmit)) and rAu(t) thus obtained are 
shown in Figs. ^ and |^, respectively. They exhibit peculiar time evolution within the 
'elongation-contraction' period. One may at once notice from r„(t) that the chain is pulled 
by the field in Intervals 2 and 3. Particularly in Interval 3, r„(t) at 40 < n < 120 becomes 
about 10~^ whose integrated weight reaches about 10%, i.e., among 10 extended parts one of 
them is as large as several tens. The latter is quite elongated and is considered to constitute 
the longer arm of the U-shaped chain. Interestingly, however, we can see from Fig. ^ that 
f/iu{t) at larger than ten is always negligibly small. This means that sliding dynamics 
of the chain even in Interval 3 consists of a large number of defect movements whose Am are 
less than several times of the mean bead distance. Thus the dynamics even in the elongated 
portion of the chain is considered to be quite stochastic when it is looked at through Am, 
i.e., changes in a unit of time of our BD simulation. In this context, we note here that up 
to several hundreds multi-scattering processes between beads are involved in one MD step 
of our BD method as explained in I. 

Another interesting observation is that both r„()f:) and rAn(^) in Interval 4 is the most 
closer, among the four intervals, to a Poisson distribution which is observed in equilibrium 
under E = Q. Namely, although Ri{t) is still relatively larger in this interval just after the 
chain is released at t tmax from trapping by an obstacle(s), many s-beads have invaded 
into the extended part(s) and the chain contracts dominantly by the entropic effect. When 
Ri{t) comes closer to its local minimum value due to trapping by a new obstacle(s), i.e., 
in Interval 1, r„(t) and rAn(^) at larger n and Au become bigger than those of Interval 4, 
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indicating that some parts of the chain are already elongated by the field. 

Now let us turn to Figs. ^ and |^ where r„ and r^u averaged over the whole time window 
of the BD run under different fields are shown. Up to the field strength of ~ 0.008 
neither r„ nor r^u significantly deviates from a Poisson distribution, while in r„ under 
E > 0.016 a tail clearly shows up at large n. This clearly indicates the existence of a 
crossover between dynamical regimes with and without the elongation-contraction motion 
at around E ^ 0.008 with M = 160 and Ogei = 20. Even in the regime with the elongation- 
contraction motion where r„ exhibits a significant deviation from a Poisson distribution, 
Am remains in an exponential distribution, though its width moderately increases with 
increasing E. This result again indicates that, even in the deterministic regime where the 
time evolution of averaged Ri{t) or VQ{t) looks apparently deterministic, defect dynamics is 
still purely stochastic. 

V. EXTENDED BOND FLUCTUATION METHOD 

The stochastic nature of 'defect' dynamics found by the BD method in the present work 
was just the basic assumption a priori introduced in the extended bond fiuctuation method 
(EBFM) when we developed it.lil The conventional BFM (CBFM)0 is a well-established 
MC method which simulate thermodynamic properties of polymers by a coarse-grained 
chain model on a lattice. It consists of a set of local movements which suffice the conditions 
necessary for the chain to evolve by self-avoiding random walk. In addition to this local 
updating process we introduced non-local movements of 'defects' in the EBFM. For this 
purpose 'defects', or 's-monomers' are defined in almost the same way as s-beads are defined 
in the present work. We then let s-monomers hop non-locally within an extended portion 
of the chain that each s-monomer lies. The movements are incorporated in such a way that 
they fulfill the detailed valance and the self-avoiding conditions. In the EBFM this non-local 
movement of s-monomers and the local updating of the CBFM are tried alternatively. We 
note that the number of s-monomers does not change except at the chain ends by the former 
process, but it does change anywhere in the chain by the latter updating. 

In Fig. ^we show r„ and of a chain with M = 80 obtained in the EBFM under = 0, 
where we denote the normalized field as = qE\^/k-BT for the MC analysis. As is the case for 
the BD model, both r„ and obtained by the EBFM obey a Poisson distribution except 
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for regions n, Au < 5. Furthermore, in the EBFM, their widths coincide with each other at 
large n and Au. In this context, it is noted that we have adopted equal probability of semi- 
local movement of an s-monomer by any distance k among possible n ones. The width of 
the distribution r„ of the EBFM is about a half of that of the BD method. Correspondingly, 
'defect' density p (~ 0.33) of the EBFM is about twice larger than that of the BD method. 
Although there exist such quantitative differences between the two methods, the qualitative 
agreement in the results demonstrated in Figs. ^ and |^ supports our EBFM approach in 
Ref. |5. 

The elongation-contraction motion under CFGE is also observed by the EBFM. The 



normalized histograms r„(t) and r/^u{t) obtained in the four intervals introduced in Sec. IV 
are shown in Figs. ^ and |10|. Their overall features qualitatively agree with those obtained 
by the BD calculation. By a detailed inspection of Fig. ^, however, we see that r„(t) in 
Intervals 2 and 3 by the EBFM still exhibit nearly exponential distributions at larger n in 
contrast to those of the BD results which have an enhanced tail at n as large as several 
tens. This difference is considered to originate from the following process of the EBFM: 
s-monomers are created or annihilated, and so length n of extended parts are changed, only 
through local updating of monomers originally incorporated by the CBFM. By such local 
stochastic processes an extremely long extended part is hardly created. 

The normalized histogram r/^u{t) shown in Fig. |rD|, on the other hand, exhibit exponential 
behavior at Au > 2 quite similar to the corresponding BD result shown in Fig. ^. This is also 
the case for the integrated distributions under various field strengths. They are shown in 
Fig. |ll] which one can compare with the corresponding BD result in Fig. 0. Quantitatively, 
fAu{t) and by the EBFM have significant weight at relatively larger Am as compared 
with those by the BD calculation. This difference is attributed to the process in the EBFM 
that, once a long extended part of length n is created, a tangential move of the s-monomer 
by the distance n is accepted by the probability proportional to 1/n. Although there exist 
such minor differences between the dynamical processes involved in the EBFM and the BD 
method, the results demonstrated above are sufficient for us to conclude that the EBFM 
properly takes into account the stochastic dynamics of slack segments which give rise to the 
elongation-contraction motion of a polymer under gel electrophoresis. 
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VI. CONCLUSION 



In equilibrium under a vanishing field it is found, as expected, that s-beads, or 'defects', 
are distributed randomly. Both distributions of extended parts r„(i) and of their tangential 
movements r/^u{t) obey a Poisson distribution. In the relatively strong field region where 
a chain exhibits the elongation-contraction motions giving rise to quasi-periodic behavior 
in Ri{t) and VQ{t), rn{t) is observed to show the corresponding time evolution. Especially, 
immediately before Ri{t) reaches to a local maximum, chain conformations with extended 
parts of more than a few tens of beads, which are entropically unfavorable, are observed. 
As for rAu{t): such a large fiuctuation seems not to occur. We therefore conclude that the 
s-bead dynamics at a semi-microscopic level is rather stochastic, though the coarse-grained 
quantities such as Ri{t) look apparently deterministic during the elongation-contraction 
motion, or in time ranges which wc have called deterministic ones. We have also shown 
that these features of r„(t), rAu{t) and R-iit) are similar to those observed by the MC 
analysis based on the EBFM which we recently proposed. Such coincidences between the two 
simulations strongly support our algorithm of EBFM which requires an order of magnitude 
less CPU times than MD calculations based on the present BD method do. 
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Figures 



FIG. 1: The normalized histograms r„ and of a chain with M = 160 under E = {) for various 
k-^T obtained by the BD calculation. 

FIG. 2: r^Au versus k^T. 

FIG. 3: The mean peak shape Ri(t — t^^^) / Riit-^iax) of — 160 with Ogd = 20 under E = 0.032. 
In each of the four intervals indicated by the arrow, histograms of Nn{t) and N/\u{t) of a sequential 
2000 steps are accumulated. 

FIG. 4: Normalized histograms r„(t) represents characteristic instantaneous distributions of n for 
the intervals defined in Fig. ^ 

FIG. 5: Normalized histograms r/^^it) in the four intervals. 
FIG. 6: r„ by the BD calculation. 
FIG. 7: TAu by the BD calculation. 

FIG. 8: The normalized histograms and r^u under = obtained by the EBFM. 
FIG. 9: The normalized histogram r„(t) obtained by the EBFM. 
FIG. 10: The normalized histogram rAM(t) obtained by the EBFM. 
FIG. 11: tam in the EBFM with various field strength O. 
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FIG. 3 
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FIG. 4 
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FIG. 10 




Au 



24 




25 



